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Using a string inspired model we construct gravity duals generalizing p x and p x + ip y supercon- 
ductors. Introducing a Chern-Simons coupling in the gravity side we demonstrate the ability to 
control which phase dominates at low temperatures, and focus on the chiral p x + ip y phase. We 
study the fermionic spectral function and establish that the behavior is rather different from the 
standard p-wave two-nodes model. 



Introduction — The discovery of materials with pos- 
sible chiral p-wave superconductivity has gathered much 
interest in condensed matter physics. This type of un- 
conventional superconductors [1-3] features a complex 
vector order parameter (abbreviated as the p + ip wave), 
time reversal breaking, and a complex gap function with 
the phase acquiring a 2w shift when rotated around the 
symmetry axis, i.e., A±(k) oc (k x ± ik y ). In theoretical 
models p+ip superconductivity is also related to electron- 
hole symmetry and the breaking of spin rotation. Chiral 
p-wave superconductivity was reported in Sr 2 Ru04 [4, 5], 
some heavy fermion compounds [6, 7], and was suggested 
in the study of of the Iron-pnictide superconductors (SC) 
[8-10]. 

The AdS/CFT correspondence [11-13] provides calcu- 
lational tools for approximations in strongly coupled field 
theories, inspiring a search for connections between string 
theory and low energy physics of compound materials like 
the cuprates, graphene, non-Fermi liquids and more. 

Review — There are three main models to study un- 
conventional superconductivity via the AdS/CFT dual- 
ity: The Abclian Higgs model [14-16], the non-Abelian 
model [17, 18], and the charged spin 2 model [19-21]. The 
distinctive feature of these models is the order parame- 
ter symmetries: s-wave, p-wave and d-wave respectively. 
Each of the models has its advantages and disadvantages 
both on the pure theoretical side and when applied to 
real condensed matter systems. These models and their 
deformations span the dictionary between string theory 
and superconductors in condensed matter physics. 

The p-wave model exhibits a vector order parameter 
and two superconductor phases, with order parameters 
(p x ) and (p x + ip y ). In the original model, and all previ- 
ous generalizations, the non-chiral (p x ) is the only ther- 
modynamically stable phase. In this letter we focus on a 
natural, from the string theory point of view, extension 
of the non-Abelian holographic superconductor which ex- 
hibits p + ip order parameter as the thermodynamically 
dominant phase. The extension is based on the introduc- 
tion of an anomalous global symmetry in the field theory 
which corresponds to a Chern-Simons term on the grav- 
ity side. This particular modification descends directly 
from considering gauged supergravities and has already 
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shown its potential as a description of a quantum critical 
point [22]. 

The (p x ) and (p x + ip y ) Holographic Supercon- 
ductors — The p-wave holographic superconductors 
studied in [17, 18, 23-26] are framed in a phenomeno- 
logical model that contains a dimensionlcss tuneable pa- 
rameter a measuring the relative strength between the 
AdS radius and the gauge group coupling. The unbroken 
phase is the Reissner-Nordstrom-Anti de Sitter (RNAdS) 
black hole. At zero temperature this is the extremal so- 
lution. To identify the potential instabilities it is worth 
discussing the AdS2 geometry arising in the near hori- 
zon region of the extremal solution. The competition 
controls the effective mass of a non-Abelian mode in 
the near horizon geometry. A strong enough YM cou- 
pling can generate a negative effective mass below the 
Breitenlohner-Freedman (BF) bound of A0IS2 and thus 
drive an instability and a global phase transition. 

Gubser's model in D-dimensions can be summarized by 
the Einstein- Yang Mills action with SU{2) gauge group, 
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with i,j = 1,2,3 and the cosmological constant A^ = 
_(Z5-i^d-2) There are two possible types of non- 
Abelian condensate distinguished by the value of the 4- 
form FL-. A FL-.. For the (p x ) case the 4-form vanishes, 

while for the {p x + ip y ) it has a finite value. The gauge 
field equation of motion prohibits turning on these two 
modes simultaneously. 

We modify Gubser's model in 5 dimensions by intro- 
ducing a Maxwell field and a Chern-Simons term: 



5* — So + S\ 
1 

8^: 



G (2 ) A *G (2) + 1 F( 2) A Ff 2) A B {1) 



G(2) =dB 



(!)• 



(3) 

1 

(4) 
(5) 



The form of the modification is inspired by 5-dimension 
SU(2) x U(l) gauged supergravity (which is a consistent 
reduction of type IIB string theory). However, the full 
supergravity contains additional neutral scalars, and our 
constants 7 and a are replaced by scalar functions. 
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On the gravity side, the Chern-Simons term induces an 
explicit breaking of the time-reversal symmetry. In the 
holographic field theory dual the Chern-Simons term in- 
troduces a 'tHooft anomaly between the SU(2) and U(l) 
global currents. The parameter 7 measures the strength 
of the anomaly and it is fixed by the dynamical properties 
of the field theory. In order to interpret these solutions as 
superconductors, one needs to weakly gauge the Cartan 
of the SU(2), i.e. (this is the field theory photon). 

The anomaly will then break the U(l) global symmetry, 
and it cannot be identified as a conserved current in the 
IR field theory. The additional two non-Abelian currents 
(-4(1) an d 4^)) remain as conserved currents in the IR. 
The appearance of the SU(2) currents in the field theory 
forces electron-hole symmetry for any charged fermionic 
excitation. Alternatively, the currents can be interpreted 
as enhancement of electron-hole symmetry in the low- 
energy effective action. 

Note that the (p x ) solution is simply embedded in the 
new model in a 7 independent manner. Since FL^ A 

FLs = there will be no source to Maxwell field 



the probe limit, the general equations become 



and without external sources it will vanish G 
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We can gain intuition into the effect of the Chern- 
Simons term by studying the model near an AdS2xR 3 
throat (the near horizon of an extremal RN solution). 
The solution Ansatz is 
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Where ((p) and u(p) are respectively the (p x ) and 
(Px + ipy) modes. The linearized equations of motion 
for p(p), oj(p) and ((p) 

p"(p)=o, 

[A'W]' + y(W=0, 
[pV(r)]' + ya;(p)=0. 
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Notice that the only parameter governing the equation 
is the dimensionless parameter a. We find an instability 
for both non-Abelian modes (comparing the mass to the 
BF bound of AdS 2 ) if 



a > a c = V3/2. 
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Remaining in the probe limit, adding the non- linear 
terms we find 



C(rV(r) = 0, ( 2 (r)u> 2 (r) = 0, 
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which is the statement that the two condensates cannon 
coexist. Focusing on the (p x + ip y ) condensate, still in 
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For non-zero 7 an instability of u>(p) will source the 
Maxwell field which in turn backreacts on oj{p). There 
is no reliable probe limit for T = 0, so to find an honest 
solution we need to take into consideration the equation 
for the metric and the charge density A^y 

The most urgent question is whether the Chern-Simons 
term can change the thermodynamics of the system. 
There is no analytical answer for this question and we are 
forced to turn our focus to numerical solutions. When 
discussing thermodynamic stability it is important to 
consider the full solution (not only the probe limit), as 
back-reaction can change the nature of the phase transi- 
tion [23] . 

Numerical solutions - We solve the equations of motion 
numerically (beyond the probe limit), and scan the effect 
of 7, for simplicity we set the AdSs scale L = 1. We use 
the following Ansatz 
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The leading terms of the boundary expansions are: 
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Uq and U\ are (proportional) to chemical potential and 
charge density of the solution. We used the scaling sym- 
metries to set the asymptotic value of G, H and / to 
1. We also set the non-normalizable modes of £, lo and 
p to zero leaving the corresponding normalizable modes 
Zi, W\ and P\ free. Finally there are two additional 
normalizable modes M\ and G\ related to metric terms. 
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The leading terms of the near horizon expansions are: geometry at a large radius 
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The parameters and St 



are related to scaling symmetries, the scaling is fixed at 
the boundary expansion so we need to follow it in the 
near horizon expansion (and cannot fix it). Regularity 
in near horizon eliminates potential log(r — r^) terms in 
the expansion. The elimination of a constant term for 
u(r) fixes the gauge choice, to the gauge where Uo is 
interpreted as the chemical potential. The vector hair at 
the horizon is given by Co, Wo and pq. The three types 
of solutions we consider are: (1) lu = p n = (2) Co — 
and (3) Abelian solution with all hair vanishing. These 
solutions correspond respectively to (p x ), {Px + ipy) and 
no condensate phases. 

We solve the equations by using a two section shoot- 
ing algorithm, matching the near horizon and the AdSs 
boundary expansions. We use numerical solutions in the 
intermediate sections to integrate the equation for a near 
horizon initial radius r, and asymptotic initial radius r/ 
to an intermediate value r m . At the intermediate point 
we define a multi-dimensional goal function as the dif- 
ference between all variables (and first derivatives). The 
goal function is a 12-dimensional function of the 14 ini- 
tial condition parameters at r.j and rf. Using a Newton- 
Raphson algorithm we find a root of the goal function 
fixing all parameters as functions of the black hole chem- 
ical potential, the black hole temperature and the La- 
grangian parameters a, and 7. We further set the chem- 
ical potential such that Uq = 1, thus working in a grand- 
canonical ensemble and measuring all quantities in units 
of the chemical potential. 

Regularization of the Euclidean action - The regular- 
ized Euclidean action is [27] 
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where the Euclidean bulk action is the analytic continu- 
ation of (3). The surface term is the Gibbons-Hawking 
boundary term. For a Yang-Mills Einstein Gravity the 
only relevant counter term is related to the induced met- 
ric on the boundary since the solution of interest has no 
divergence coming from the Yang-Mills part. 
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with h the induced metric on the boundary and n is an 
outward pointing unit normal vector to dM. Cutting the 
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Results - We find numerically that G\ = P\ = 0, so 
the non-vanishing normalizable modes are related to the 
energy density(Ml), charge density([/i) and the two pos- 
sible condensates Z\ and U\. Figure 1 presents a typical 
thermal behavior of the solutions. We find that increas- 
ing the value of 7 lowers the free energy of the (p x + ip y ) 
solutions until the point where they win the thermody- 
namic battle against the (p x ) solutions. We also see that 
large values of 7 (around 2.5) can change the nature of 
the phase transition from second order to first order. We 
note that at smaller values of a (than those shown in 
plots) the (p x ) solution also exhibits first order phase 
transition, while the (p x + ip y ) solution is always second 
order at 7 = 0. 

The summary of the dependence of the phase diagram 
on 7 and a is given in 2. On the right side we see the 
dependence of the critical temperature on a. The left 
side holds the answer for the most urgent question, we 
see that for any value of a (in the range tested) there 
exist a critical value of 7 above which the (p x + ip y ) is 
the thermodynamically favored phase. We define j c by 
comparing the free energy at T — 0.5T C , using T c of the 
(p x ) solution. The error bars on the plot are the result 
of the heavy computation time involved in finding each 
value, and not machine precision. 

The near horizon Ansatz we use is valid for finite tem- 
perature solutions only, leading to a minimal tempera- 
ture below which we loose numerical accuracy. Therefore 
the question of which phase dominates at T = remains 
open. However for large enough values of 7 it is likely 
that the (p x + ip y ) phase will remain thermodynamically 
favored all the way down to zero temperature. 
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FIG. 1. Dependence of the condensate and free energy on 
temperature (arbitrary units). Continuous green line is the 
regular phase, dotted red line is the (p x ) phase, and dashed 
blue lines are {p x + ipy)- Plots are for a = 1.6 and the dashed 
curved realetd to 7 = 0, 7 = 1.25 and 7 = 2.5. 

The Spectral Function — Angle Resolve Photoemis- 
sion Spectroscopy (ARPES) allows to directly probe the 
fermionic structure of materials. In our theoretical model 
the fermion spectral function can be read from the re- 



FIG. 2. On the right side we plot the dependence of the 
critical temperature T c on the gauge coupling strength a. The 
continuous red line is the (p x ) phase, and the dashed blue lines 
are (p x + ip y ) with 7 = 0, 2, 2.75, 3.5. On the left we plot the 
dependence of 7 C on a. The critical value 7 C is defined as 
the minimal value of 7 where the (p x + ip y ) dominates at 
T = 0.5T C (we use the critical temperature of {p x })- 
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where Li, k are frequency and momentum in the dual field 
theory. 

We use a "minimally coupled" fermion with zero mass 
transforming in the 2 of the SU(2) gauge symmetry and 
with charge q = 1/2 under the U(l) gauge symmetry. 
See appendix-A for details. Rotations around the axis of 
symmetry (z) are canceled by SU(2) rotations, and after 
the trace the spectral function is invariant in the x — y 
plane. We use polar coordinates in the x — z plane k z — 
kcos9. We calculated numerically the spectral function 
for a = 1.59, 7 = 2.5 and temperature T/T c = 0.2. The 
results are displayed in figures 3-5. For comparison in 
appendix-B we reproduce the same results for the (p x ) 
SCs. 

In figure-3 we can see the general form of the spectral 
function in the k — uj plane at a fixed angle (k x /k z = 
tan 8) . The spectral function is invariant under reflec- 
tion in the frequency oj — > —uj. Note that this symmetry 
is merely a consequence of the trace but in general it will 
be broken. We can clearly see two separate excitations 
(resonances) , each appearing twice (negative and positive 
frequency). The Fermi-momentum (fcp) for each excita- 
tion is the momentum where the excitation frequency is 
closest to zero. Half the frequency distance between the 
positive and negative excitation (of the same type) at the 
Fermi-momenta is the excitation gap [28] . 

Figure-4 follows the change of the excitation gap and 
Fermi momentum for both excitations as we change the 
angle 8. The higher frequency excitation is gapped at all 
angles, and the gap size changes mildly with angle, while 
the corresponding Fermi momentum (k F t9h ) moves from 
zero near 8 = to 0.6 at 8 = it. The lower frequency 
excitation is not gapped at 8 = and develops a gap as 
the angle increases, the corresponding Fermi momentum 
(kp w ) changes mildly with 8. As can be seen from left 
plot in figure-3 the lower frequency excitation, can be 
strongly attenuated at k l p w . We do not have a clear 



FIG. 3. A density map of the spectral function at fixed direc- 
tion in the x — z plane. The color scheme changes from near 
zero values in blue to high values in red. The plot is clipped 
below the maximal value of the spectral function to allow to 
see the details better. On the left we look at 9 = ^ and 
6 — 5 on the right. 



explanation for this behavior, and it may be the result of 
the interplay between the two resonances. 
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FIG. 4. The excitation gap (left) and Fermi momentum 
(right) for the two excitations as a function of angle 9. The 
blue circle are related to the lower frequency excitation and 
the red triangles for the high frequency excitation. The lower 
frequency excitation (blue circles), shows a strong attenuation 
near it's Uf for a wide range of angles. When the attenuation 
is below 10 -3 (arbitrary cutoff) we replace the circles by dots. 

To further study the angle dependance of the spec- 
tral function we study the Energy Distribution Curve 
(EDC), which is defined as the spectral function at the 
Fermi-momentum. In Figure-5 we plot the EDCs for 
both Fermi-momenta. The plot on the left, which is the 
EDC at k l p w , shows that the peak is strongly attenuated 
for angles near the poles (8 = and ir). The EDC at 
k'p 9h show the high frequency peak disappears at mid 
range angles, however this is not due to attenuation near 
the Fermi-momenta, but due to the dominace of the lower 
frequency peak at these range of 9. 

In appendix-C we elaborate upon additional properties 
of the spectral function. First on the isospin current and 
then on the node structure. The isospin current is just 
the consequence of the relation between spatial rotation 
and rotation of the SU(2) group. Rotation of the SC so- 
lution in the x-y plane is not an exact symmetry, but can 
be compensated by a global rotation in the gauge group 
SU(2). This property can be detected in the fermionic 
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FIG. 5. The EDC (spectra function at the Fermi-momenta) 
as function of frequency. The different plots are related to 
different angles 8. The left and right plots are the EDC at the 
k l p w and k F ' 3h correspondingly. The dashed red line represent 
the location of the gap. The different EDO's are normalized 
differently between the angles (for this plot only) and should 
not be directly compared. 



Green's function by considering 
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The node structure is studied by examining the low fre- 
quency spectral function, for the (p x + ip y ) we see a sin- 
gle node near k — (0,0, —0.57), 

Discussion — To summarize, in the first part we 
showed how a holographic chiral p + ip superconductor 
can be constructed using the full treatment of gravity 
equations of motion. In the second part we study the 
fermionic spectral function which diplay some unusual 
behavior compared to the p-wave superconductors (see 
appendix-B), the new features found: 

• Two competing modes versus a single dominant 
mode in the (p x ) SCs 

• Chirality, i.e relation between spatial and SU(2) ro- 
tations which is absent in the (p x ) SCs. 

• Single node spectral function. 

It would be interesting to pursue some of the follow- 
ing directions. For example, it would be desirable for 
condensed matter applications to consider similar con- 
structions where the time-reversal symmetry is broken 
spontaneously. However, this would likely require a non- 
radial solution on the gravity side. 

Recently, Son and Surowka [29] have analyzed the hy- 
drodynamics in the presence of an anomalous current and 
its description in terms of gravity. It would be interest- 
ing to use their approach to analyze our gravity solutions 



which should lead to the hydrodynamics of a p + ip su- 
perconductor. 
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Appendix A: Calculation of the spectral function 

The tree level contribution for the Green function 
is calculated by solving the Dirac equation of a probe 
fermion in the black hole background with in-falling 
boundary condition at the horizon. We use a "minimally 
coupled" fermion with mass M transforming in the 2 of 
the SU(2) gauge symmetry and with charge q under the 
U(l) gauge symmetry. 
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where T a are flat space Dirac matrices, E£'s are the in- 
verse Vielbein, and r"s are generators of the su(2) alge- 
bra. We use the following redefinition of the fermion field 
= (-gg rr )~^e- iuot+ik i xi ipi. Near the AdS 5 bound- 
ary the Dirac equation can be solved in terms of 4 nor- 
malizable and 4 non-normalizable modes 
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where Aj a , Dj^ are free boundary conditions and Bj a , 
Cia are constants fixed from the equation of motion. 
Note that v a and u a are eigenvectors of T- with eigen- 
value +1 and —1 respectively. 

In the near horizon limit (finite temperature) the met- 
ric takes the form 



ds 2 = -f?(r-r h )dt 2 + 



Prdr 2 



dx 2 
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The Dirac equation solution behaves like 

4 

1>i = J2 [ a mVr{r - r h y^ + b m U'j n (r - rrf^ 

m— 1 

(A4) 

where Vf 1 and U™ are a vector base for in-falling and 
out-going solutions respectively and a m and b m are free 
parameters. Since the equations are first order, setting 
in-falling boundary condition (b m = 0) impose a rela- 
tion between the asymptotic boundary conditions. The 
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asymptotic values are then determined (with a slight 
abuse of notations) as 

A = Aa, D = Da (A5) 

The retarded Green function is the 4x4 matrix. 

G R = — zD A _1 r* (A6) 

Notice that T* should be projected onto the spin states, 
and it should not multiply the SU{2) part of the Green 
function. 
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Appendix B: Spectral function of the (p x ) 
superconductors 

In this appendix we reproduce the result for the spec- 
tral function of the (p x ) SCs in order to allow for a clear 
comparison between the two types of p-wave supercon- 
ductors . To the best of our knowledge previous work 
[30] covered only the probe limit, and the presentation of 
the result is different and harder to compare to the result 
of this paper. 

In figure 6 we display the spectral function for a (p x ) 
superconductor at various directions in the momentum 
plane. Figure 6 shows four quasi-particle excitations; 
two with positive frequencies and two with negative fre- 
quencies. Following the peak in the spectral functions 
we can find their dispersion relations tu(k). We can esti- 
mate their individual contributions to the whole spectral 
function by integrating the spectral function along these 
curves (spectral weight of a single excitation). There 
are also nodes near 9 = 0, it where there is no gaps in 
the spectral function J k p(0, k) > 0. One can see the 
the inner excitation (closer to k = 0) is dominating for 
all the angles, in the sense that it has higher spectral 
weight. Also one can see clear gapless excitation near 
the pole 9 = 0,ir. The spectral function show symmetry 
under rotation in the y — z plane, and flipping of the fre- 
quency (a; — > — lu) and momentum along the symmetry 
axis (k x — > —k x ). For comparison we plot the matchings 
plots of the (p x + ip y ) superconductor in figure 7. The 
qualitative difference is very clear, most prominently is 
the existence of two clear nodes near 9 = and it. A 
second clear difference is invariance of the (p x ) spectral 
function under 9— > n — 9 (flipping the sign of the mo- 
mentum). A more subtle difference is the dominance of 
the lower frequency excitation in the (p x ) case, in con- 
trast to the (p x + ip y ) spectral function where the two 
excitation alternate dominance with the angle and mo- 
menta (most obvious is the nearly disappearance of the 
lower frequency near 6 = vr). 



Appendix C: Additional aspects of the Green's 
function 



In this section we give extra (numerical) details on the 
nature of the fermionic excitation as can be seen from 



FIG. 6. The spectral function of a (p x ) superconductor as a 
function of momentum (horizontal axis) and frequency (verti- 
cal axis). Plots corresponds to different angles 6 = tan" 



-1 £x 



The angle used on the top line from left to right are are 
= 0, y , and . The angles used on the bottom line 
from left to right are 9 = ^f, ^f, ^ and tt. Unlike figure-3 
this plot use a logarithmic scale for the colors (see color table 
below the plot), which help to better compare the height of 
different peaks. The momentum range used is from to 1 the 
frequency range used is from — 1 to 1 . 
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FIG. 7. The spectral function of a {p x 
See the caption of figure 6 for details. 



■ ip y ) superconductor. 



the 2-point function. 

The fermionic 2-point function is a 4 x 4 matrix, where 
the 4 is combination of the global SU(2) symmetry and 
(/ = 1,2) and the projection over r r eigenvalues (a = 
1,2). So far we studied the spectra function which is the 
imaginary par of the trace: 
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The large SU(2) and spatial rotations transformations are 
non-gauged and appear as global symmetries in the field 
theory. The conserved currents corresponding to these 
symmetries is measurable in Josephson junctions exper- 
iments [31, 32]. Out of the possible currents only two 
are not vanishing, the isospin and spin current densities. 
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They are defined as follows: 
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FIG. 8. The isospin (left) and spin (right) currents. The 
blue circles mark the {p x + ip y ) superconductor, while the 
red squared mark the (p x ) superconductor. The angle <f> is 



the angle around the symmetry axes tan 



and tan 



for the {p x + ip y ) and (p x ) superconductors respectively. The 
plot are taken at momentum | k | =0.88. 



/(<**»)(£) = - - J2 [Tr Gfl(w, k)iT 2 



la, 1/3 



The local density of states I(k) = J^u, p{ u i w) is invari- 
ant under rotations in the symmetry plane. However the 
isospin and spin currents densities change as we rotate 
around the symmetry axis. In figure-8 we plot the angle 
dependence of these currents for the {p x ) and (p x + ip y ) 
superconductor. As expected from the symmetries of 
the background, the isospin current vanishes for the (p x ) 
type, but show a perfect sin behavior for the (p x + ip y ) 
SC. The spin current is non-trivial for both cases. 




Another interesting view, that give a better description 
on the location of the nodes is the spectral function in the 
momentum plane for fixed frequency. See figures and for 
the two type of the superconductors. For the {p x ) we can 
see how the two nodes (peaks), appearing at the poles for 
low frequency. As the frequency increase the peaks move 
toward the center, and the spectral function turns into 
a Dirac cone. For the (p x + ip y ) superconductor we see 
that the low frequency picture is a single node at 9 = tt. 
As the frequency increase the node split into three peaks 
before developing into a Dirac cone. 
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FIG. 9. The spectral function of a (p x ) superconductor 
as function of momenta for fixed frequency. Going from 
left to right the top row plots are for frequencies are ui = 
0.025,0.05,0.075,0.01,0.125, and the bottom row plots fre- 
quencies are u) = 0.15,0.225,0.3,0.375 and 0.45. We use a 
logarithmic color scale as described on the bottom of the plot. 
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FIG. 10. Spectral function of a {p x + ip v ) SC as function of 
momenta for fixed frequency. See the caption of figure-9 for 
details. 



